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ABSTRACT 

A three-dimensional, Newtonian hydrodynamic technique is used to follow the postbounce phase of a 
stellar core collapse event. For realistic initial data we have employed post core-bounce snapshots of the 
iron core of a 20 Mq star. The models exhibit strong differential rotation but have centrally condensed 
density stratifications. We demonstrate for the first time that such postbounce cores are subject to a 
so-called low-r/|IU| nonaxisymmetric instability and, in particular, can become dynamically unstable 
to an TO = 1 - dominated spiral mode at T/|IU| 0.08. We calculate the gravitational wave (GW) 
emission by the instability and find that the emitted waves may be detectable by current and future 
GW observatories from anywhere in the Milky Way. 

Subject headings: hydrodynamics - instabilities - gravitational waves - stars: neutron - stars: rotation 


1. INTRODUCTION 

Rotational instabilities are potentially important in the 
evolution of newly-formed proto-neutron stars (proto- 
NSs). In particular, immediately following the pre¬ 
supernova collapse - and accompanying rapid spin up - 
of the iron core of a massive star, nonaxisymmetric insta¬ 
bilities may be effective at redistributing angular momen¬ 
tum within the core. By transferring angular momentum 
out of the centermost region of the core, nonaxisymmet¬ 
ric instabilities could help explain why the spin periods of 
newly formed pulsars are longer than what one would ex¬ 
pect from standard stellar evolutionary calculations that 
do not invoke magnetic field action for angular momen¬ 
tum re distribution and generation of very slowly rotating 
cores llHeger et al.l l200(ilHirschi et all l2004lHeger et al.! 
1200411 . Alternatively, in situations where the initial col¬ 
lapse “fizzles” and the proto-NS is hung up by centrifugal 
forces in a configuration below nuclear density, a rapid 
redistribution of angular momentum would facilitate the 
final collapse to NS densities. The time-varying mass 
multipole moments resulting from nonaxisymmetric in¬ 
stabilities in proto-NSs may also produce GW signals 
that are detectable by the burgeoning, international ar¬ 
ray of GW interferometers. The analysis of such signals 
would provide us with the unprecedented ability to di¬ 
rectly monitor the formation of NSs and, perhaps, black 
holes. 

In this Letter, we present results from numerical sim¬ 
ulations that show the spontaneous development of a 
spiral-shaped instability during the postbounce phase 
of the evolution of a newly formed proto-NS. These 
are the most realistic such simulations performed, to 
date, because the pre-collapse iron core has been drawn 
from the central region of a realistically evolved 20 
Mq star, and the collapse of the core as well as the 
postbounce evolution has been modeled in a dynami¬ 
cally self-consistent manner. Starting from somewhat 
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simpler initial states, other groups have followed the 
development of b ar-like structure in p ostbounce cores 
using Newtonian lIRamp p et alJ 1199811 and relativistic 
(IShibata & Sekiguchil I2005ll gravity, but their analyses 
have been limited to cores having a high ratio of rota¬ 
tional to gravitational potential energy, f3 = T/|IU| > 
0.27. We demonstrate that a one-armed spiral (not the 
traditional bar-like) instability can develop in a proto- 
NS even if it has a relatively low r/|IU| ~ 0.08. This is 
significant, but perhaps not surp rising given the recent 
studies by Ceni/enaoet ajj j200 ill . Shibata et al. (2002, 
2003), and Saiio et al.1 ll2003tl . 

2. NUMERICAL SIMULATION 

The results presented in this Letter are drawn from 
three-dimensional hydrodynamic simulations that follow 
approximately 130 ms of the “postbounce” evolution of a 
newly forming proto-NS. Before presenting the details of 
these simulations, however, it is important to emphasize 
the broader evolutionary context within which they have 
been conducted and, specifically, from what source(s) the 
initial conditions for the simulations have been drawn. 
The two simulations presented here cover the final por¬ 
tion (Stage 3) of a much longer, three-part evolution that 
also included: (Stage 1) the main-sequence and post- 
main-sequence evolution of a spherically symmetric, 20 
Mq star through the formation of an iron core that is 
dynamically unstable toward collapse; and (Stage 2) the 
axisymmetric collapse of this unstable iron core through 
the evolutionary phase at which “bounces” at nuclear 
densities. 

Stage I of the complete evolution was originally pre¬ 
sented as model “S'20” by Woosley & Weaver (1995). 
The initial configuration for this model was a chemically 
homogeneous, spherically symmetric, zero-age main- 
sequence star with solar metallicities. Evolution up to 
the development of an unstable iron core took some 
2 X 10^ yr of physical time. In Stage 2 the spheri¬ 
cally symmetric model from Stage 1 was mapped onto 
the two-dimensional, axisymmetric grid of the hydrody¬ 
namics code “VULCAN/2D” (Livne 1993) and evolved 
as model “5'20A500/30.2” by Ott et al. (2004). Rotation 
was introduced into the core with a radial angular ve¬ 
locity profile = Do[l -I- (w/AY]~^ (where w is the 
cylindrical radius). The scale length in the initial rota- 
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tion law was set to ^ = 500 km, and flo = 3.36 rads“^ 
was chosen so that, initially, [3 = 0.0020. The axisym- 
metric collapse was modelled adiabatically but the full 
Lattimer-Swesty equation of state (LSEOS; Lattimer & 
Swesty 1991) was incorporated. 

During Stage 2, the innermost region of the unsta¬ 
ble iron core collapsed homologously and, in « 0.5 s, 
reached nuclear densities. As a consequence of angular 
momentum conservation, the core spun up considerably; 
at the time of the bounce, the rotational energy param¬ 
eter had increased to [3b = 0.0896. As a result, the core 
bounce was aided as much by increased centrifugal forces 
as it was by the rapidly stiffening equation of state at 
nuclear densities. After bounce, the core expanded co¬ 
herently, leading to almost an order-of-magnitude drop 
in the maximum density. This expansion was then re¬ 
versed when gravitational forces again began to dominate 
over pressure gradients and centrifugal forces. In this 
way, the rapidly spinning, postbounce core underwent 
several damped-harmonic-oscillator like cycles. Ott et 
al. (2004) followed the axisymmetric evolution of model 
S'20A500/30.2 for ~ 140 ms after the time of its initial 
bounce, tf,- 

In order to study the possible development of nonax- 
isymmetric structure in a newly forming proto-NS, we 
mapped the VULCAN/2D model into the fully three- 
dime nsional hydrodynamic code FLOW*ER lIMotl et al.l 
l2002t Stage 3). In this work, we have adopted an ideal- 
fluid EOS: p = (r — l)ep, where p is the mass density, p 
is the pressure, T is the chosen adiabatic exponent, and 
e is the specific internal energy. For sub-nuclear density 
matter — specifically, for p < pnuc = 2 x 10^“^ g cm“^ 
— r is set to 1.325, which approximates the T given 
by the LSEOS in VULCAN/2D in the sub-nuclear den¬ 
sity regime and the given conditions. For p > pnuc, we 
have set T = 2.0 to mimic the effects of nuclear repulsive 
forces. FLOW*ER’s uniform cylindrical grid was con¬ 
structed in such a way that it enclosed the innermost 140 
km of model S'20A500/30.2, containing around 1.4 Mq (~ 
75% of the mass that was on the VULCAN/2D grid at 
the postbounce stage). Model 520^500/30.2 was mapped 
onto FLOW*ER’s grid at two different times during its 
postbounce evolution: Model “Q15” was evolved on a 
grid with (98,128,194) zones in (tu, (/>, Z) and mapped at 
t-tb = 15 ms when the core was in the middle of its second 
postbounce expansion phase. Upon introduction into the 
three-dimensional code, its density field was perturbed 
with an 0.1% amplitude, bar-like m = 2 seed. Model 
“W5” was evolved on a higher resolution (130, 256,194) 
grid, beginning at t-tf, = 5 ms when the model reached its 
first density minimum after bounce. Random perturba¬ 
tions of 0.02% amplitude were imposed on the densities 
at mapping. Both models were evolved up to a time 
t — tb ~ 130 ms. 

3. RESULTS 

The key results of our three-dimensional evolutions 
are displayed in Fig. Q] (for model Q15) and Fig. [21 (for 
model W5). In the top panel of each figure, the time- 
dependent behavior of the global dipole (m = 1; solid 
curve) and quadr upole (m = 2; dotted curve) moments 
ijSaiio et abll200^ are plotted to illustrate how, and on 
what timescale, nonaxisymmetric structure developed. 
The time-dependent behavior of the rotational energy 



(rns) 

Fig. 1. — Time-evolution of various physical quantities is shown 
for model Q15; time (in ms) is given relative to Top: Globally- 
averaged amplitude of m = 1 (solid curve) and m = 2 (dotted 
curve) distortions. Middle: The rotational energy ratio (3 (solid 
curve) and the core’s maximum density pmax (dotted curve). Bot¬ 
tom: Product of the GW strain and the distance to the source 
r as viewed down the rotation [z) axis (solid curve) and as viewed 
along the fc-axis (dotted curve). 
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Fig. 2.— Same as Fig. 1, but for model W5. 


parameter, p (solid curve), and the core’s maximum den- 
sity, Pmax (dotted curve), are shown in the middle section 
of the figures. The bottom third of each figure displays 
the amplitude of the gravitational radiation that would 
be emitted from each model as estimated by the post- 
Newtonian quadrupole formalism (see, e.g., Misner et al. 
1973); specifically, for each model the time-dependent be¬ 
havior of the product of the “plus polarization” of the 
GW strain and the source distance r is shown as seen by 
an observer looking down the rotation axis (solid curve) 
or perpendicular to that axis (dotted curve). For a core¬ 
collapse event at r = 10 kpc, an amplitude rh^ = 100 cm 
translates into h+ = 3.2 x 10“^^. 

Because nonaxisymmetric perturbations were initially 
introduced into both models at a very low amplitude, 
the early phase of both evolutions resembled the axisym¬ 
metric evolution reported in Ott et al. (2004). Here, 
this is illustrated best by the oscillations in /3(<) and 
Pmax{t) that are shown in Fig.|21for model W5; the char¬ 
acteristic (dynamical) time between successive “radial” 
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bounces is 2r<iyn ~ 4 ms (consistent with a mean core 
density p ~ = 1.2 x 10^^ g cm“^) and the 

axisymmetric oscillations persist for ~ 50ms. The curves 
of h+{t) for model W5 also signal that the dynamics is 
essentially axisymmetric: as viewed along the cc-axis, the 
GW strain exhibits oscillations of diminishing amplitude, 
as reported in Ott et al. (2004), but for the first ^ 40 ms 
after 4, essentially no GW radiation is emitted along 
the z-axis. In model Q15, fewer “radial” bounces occur, 
they damp out somewhat more rapidly, and the resulting 
h+(t) signal is weaker as viewed along the cc-axis. This is, 
in part, because the postbounce core configuration was 
introduced into FLOW*ER at a later time (t—tb = 15 ms 
for model Q15 instead oi t — tb = 5 ms for model W5) 
and, in part, because the effects of numerical damping 
are inevitably more apparent when a simulation is run 
on a grid having lower spatial resolution. As is illus¬ 
trated by the solid h+{t) curves in the bottom panels of 
Figs. Hand El at early times the amplitude of the gravi¬ 
tational radiation that would be emitted along the z-axis 
is larger in model Q15 than in model W5. This reflects 
the fact that the nonaxisymmetric perturbation that was 
initially introduced into model Q15 was larger and it had 
an entirely m = 2 character. 

Although in model Q15 the postbounce core was sub¬ 
jected to a pure, m = 2 bar-mode perturbation when it 
was mapped onto the FLOW*ER grid, the amplitude of 
the model’s mass-quadrupole distortion did not grow per¬ 
ceptibly during the first 100 ms ('^ 50 dynamical times) 
of the model’s evolution (Fig. EJ. However, as the solid 
curve in the same figure panel shows, the model spon¬ 
taneously developed an to = 1 “dipole” distortion even 
though the initial density perturbation did not contain 
any to = 1 contribution. As early as t — tb ~ 70 ms, a 
globally coherent to = 1 mode appeared out of the noise 
and grew exponentially on a timescale Tgrow ~ 5 ms. At 
t — tb ~ 100 ms, the amplitude of this to = 1 distor¬ 
tion surpassed the amplitude of the languishing to = 2 
structure and, shortly thereafter, it became nonlinear. 
Mt — tb w 100 ms, the quadrupole distortion also began 
to amplify, but it appears to have only been following 
the exponential development of the to = 1 mode. An 
analysis of the oscillation frequency of both modes re¬ 
veals them to be harmonics of one another. As the top 
panel of Fig. El illustrates, the same to = 1 mode devel¬ 
oped spontaneously out of the 0.02% amplitude, random 
perturbation that was introduced into model W5. The 
mode reached a nonlinear amplitude somewhat earlier in 
model W5 than in model Q15, presumably because the 
initially imposed random perturbation included a finite- 
size contribution to an to = 1 distortion whereas the den¬ 
sity perturbation introduced into model Q15 contained 
no TO = 1 component. The growth timescale of the in¬ 
stability is Tgrow ~ 4.8 ms for model W5. Although we 
have described the unstable to = 1 mode as a “dipole” 
mass distortion, this is somewhat misleading because in 
neither model did the lopsided mass distribution produce 
a shift in the location of the center of mass of the system. 
Instead, as is illustrated in Fig. El the mode developed 
as a tightly wound, one-armed spiral, very similar to the 
TO = 1 - dominated stru cture s that have been reported 
bv iGentrella et all ll2001ll and lSaiio et al.l ll2003ll . 

After the spiral pattern reached its maximum ampli- 



Fig. 3.— The equatorial-plane structure of model W5 is shown 
at time t — = 71 ms. Left: Two-dimensional isodensity contours 

with velocity vectors superposed; contour levels are (from the in¬ 
nermost, outward) p/pmax =0.15,0.01,0.001,0.0001. Right: Spiral 
character of the m = 1 distortion as determined by a Fourier anal¬ 
ysis of the density distribution; specifically, the phase angle 
of the m = 1 Fourier mode is drawn as a function of 'uj. 



Fig. 4.— Equatorial-plane profiles of the azimuthally averaged 
angular velocity ^('Cij) (top frame) and the mass density p(rt7) (bot¬ 
tom frame) are shown at four different times during the evolution 
of model W5. Changes in these profiles at late times illustrate the 
effects of angular momentum redistribution by the m = 1 spiral 
mode: angular momentum migrates radially outward while mass 
migrates radially inward. A horizontal line drawn in the top panel 

^CR = 2.5 X 10^ rad s~^ identifies the corotation radius for this 
one-armed spiral mode. The “kink” seen in p(^) at late times at 
about 8 km is connected to the discontinous switch of the EOS F 

pnuc* 

tude in both of our model evolutions, the maximum den¬ 
sity began to slowly incre ase and 0 s t arted decreasing 
(Figs. 1 & 2). Following iSaiio et al.l ll2003ll . we inter¬ 
pret this behavior as resulting from angular momentum 
redistribution that is driven by the spiral-like deforma¬ 
tion and by gravitational torques associated with it. As 
angular momentum is transported outward, the centrifu¬ 
gal support of the innermost region is reduced, a larger 
fraction of the core’s mass is compressed to nuclear den¬ 
sities and, in turn, [3 decreases because the magnitude 
of the gravitational potential energy correspondingly in¬ 
creases. Fig. E] supports this interpretation. As the 
proto-NS evolves, we see that the outermost layers spin 
faster and the innermost region becomes denser. (We 
note that throughout the evolution our models conserved 
total angular momentum to within a few parts in 10^.) 
Also, as is shown in the top panel of Fig. ^ through¬ 
out most of the model’s evolution there is a radius in¬ 
side the proto-NS [wck ~ 12 km) at which the angular 
velocity of the fluid matches the angular eigenfrequency 
(wcR = 2.5x10^ rad s“^) of the spiral mode. Hence, it is 
entirely reasonable to expect that resonances associated 
with this “corotation” region are able to effect a redis- 
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tri bution of angular m om entum in the manne r described 
bv iContoDOuloa lUQSCli or lWatts et alJ l|2005D . 

Finally, we note that an off-center density maximum 
shows up at intermediate times, but since it is within 
the innermost radial grid zones, we are not sure if it is 
physical or an artifact due to the boundary conditions at 
the axis. 

4. SUMMARY AND DISCUSSION 

Using the 3D, Newtonian, hydrodynamic code 
FLOW*ER, we have modelled the postbounce phase 
of the evolution of an iron core from an evolved, 
20 M 0 Star. In our first simulation (model Q15), the 
rotating core was mapped from two dimensions onto the 
three-dimensional grid at t — tb = 15 ms and the core’s 
axisymmetric structure was altered by the introduction 
of a 0.1% amplitude, barlike density perturbation. The 
core was found to be dynamically stable to this pure 
m = 2 perturbation, but it proved to be dynamically un¬ 
stable toward the spontaneous development of a tightly 
wound, TO = 1 spiral mode. In an effort to examine the 
robustness of this result, we performed a second simula¬ 
tion (model W5) in which the rotating core was mapped 
at an earlier time (t — tb = 5 ms), onto a grid with higher 
spatial resolution, and the core’s axisymmetric structure 
was altered by adding a lower-amplitude ( 0 . 02 %) and 
random density perturbation. The to = 1 spiral mode 
developed spontaneously in this model as well. Hence, we 
conclude that even relatively slowly rotating proto-NSs 
can be susceptible to the development of a spiral-shaped, 
TO = I - dominated mass distortion. 

Although the nonaxisymmetric distortions that de¬ 
veloped in both our models did not grow to particu¬ 
larly large nonlinear amplitudes, they produced maxi¬ 
mum GW amplitudes comparable to the “burst” signal 
produced by the preceding, axisymmetric core collapse: 
Model Q15 showed the highest amplitude, ~ 

100 cm at a frequency f ~ 800 Hz. The peak amplitude 
of the axisymmetric bounce signal reported in Ott et al. 
(2004) was rh+^max ~ 300 cm at a frequency f ~ 400 
Hz. If the source is located within the Milky Way, both 
“burst” signals may be detected by the currently oper¬ 
ative GW observatories (Ott et al. 2004). In our sim¬ 
ulations, approximately 100 ms separated the peaks of 
these two GW “bursts,” but the earlier peak near t = tb 
would, in practice, be unobservable if the rotation axis 
of the proto-NS were oriented along our line of sight, as 
is likely to be the case for core collapse events associated 
with gamma-ray bursts. 

Our results demonstrate that a realistic, non¬ 
equilibrium postbounce stellar core can become dynam¬ 
ically unstable to an to = I spiral instability at a value 
of (3 as low as 0.08. (In model W5, this corre- 
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spends to a spin period of 15 milliseconds at the sur¬ 
face of the proto-NS.) The value obtained in previous 
studies llCentrella et al.ll200lL iSaiio et ^1200311 for dif¬ 
ferentially rotating equilibrium models with a consider- 
ably simpler ther modynamic structure was (3 ~ 0.14. 
ISaiio et ^ l|2003ll argue that the growth of the insta¬ 
bility in their models requires a soft EOS with effectiye 
F < 1.4. Our models exhibit an effective F just above 
this threshold while going unstable at lower (3. This dif¬ 
ference is most probably caused by the more complicated 
thermodynamic and rotational structure of our models. 
Most significantly, the instability in our low-T /\W\ mod¬ 
els appears to be related to the proto-NS’s angular ve¬ 
locity profile 13 (vu). Therefore, unlike the classical bar¬ 
mode instability that becomes unstable above a critical 
value of T/1 IF I, the relative stability of this spiral mode 
seems likely to depend on the existence or absence of a 
corotation resonance within the star. A wider variety of 
model simulations will be required to properly examine 
the validity of thi s conjecture . Our simulations support 
the suggestion of ISaiio et al.l ll2003ll that a spiral-mode 
instability can be effective at redistributing angular mo¬ 
mentum within a proto-NS. In so doing, the instability 
can spin down the bulk of the core and, simultaneously, 
assist contraction of the innermost regions toward higher 
mean densities. 

The present study marks only one very early step in 
our understanding of rotational instabilities in proto- 
NSs. The models we have used were purely hydrody¬ 
namic and Newtonian and did not include neutrino pro¬ 
duction and radiative transfer, or any of a variety of other 
microphysics that is likely to be relevant to these astro- 
physical systems. Future, fully consistent simulations in¬ 
cluding all the relevant physics will be needed to provide 
more definitive answers to these questions of stability. 
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